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Abstract. The origin of ultra-high energy cosmic rays is discussed in light of the latest 
observational results from the Pierre Auger Observatory, highlighting potential astrophysical 
sources such as active galactic nuclei, gamma-ray bursts, and clusters of galaxies. Key issues 
include their energy budget, the acceleration and escape of protons and nuclei, and their 
propagation in extragalactic radiation and magnetic fields. We briefly address the prospects for 
Telescope Array and future facilities such as JEM-EUSO, and also emphasize the importance of 
multi-messenger X-ray and gamma-ray signatures in addition to neutrinos as diagnostic tools 
for source identification. 



1. Introduction 

Several decades after their discovery, the origin of ultra-high energy cosmic rays (UHECRs), 
cosmic particles with energies 10^^-10^^ eV and above, continue to be one of the biggest mysteries 
in astrophysics [T]. During the last few years, the field has witnessed revolutionary advances, 
thanks to the advent of new generation, hybrid detector facilities spearheaded by the Pierre 
Auger Observatory in the southern hemisphere The Telescope Array |H] has begun operation 
in the northern hemisphere and we eagerly await their data. Expectations are already high 
for future projects such as Auger North [4j, or the Extreme Universe Space Observatory on 
the Japanese Experiment Module (JEM-EUSO), to be deployed on the International Space 
Station [5]. This article, by no means a thorough review, discusses selected topics in the rapidly 
developing field of UHECRs, focusing on the theory of potential astrophysical sources confronted 
with the latest observational results. Some of the issues below are covered in more detail in earlier 
articles [6]. 



2. General issues on propagation and recent observational results 

We first touch upon some general issues concerning the propagation of UHECRs (see e.g. [TJ [8] 
for more details). If UHECRs are protons of extragalactic origin, photopion interactions with 
cosmic microwave background (CMB) photons must induce severe energy losses at 6 x 10^^ eV 
for propagation lengths ^ 100 Mpc. Unless the sources lie much nearer, a spectral steepening 
(Greisen-Zatsepin-Kuzmin or GZK cutoff [9]) is expected above these energies. The Pierre Auger 
Observatory [TO] as well as the HiRes experiment [TT] have recently demonstrated the existence 
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of a strong spectral steepening at just this energy with high statistical significance, although 
at present it cannot be excluded that this represents instead some limit to acceleration at the 
sources. 

The absence of super-GZK events reduces the motivation for the so-called top-down scenarios, 
in which UHECRs arise as decay products of higher energy entities envisaged in non-standard 
particle physics models, such as topological defects, superheavy dark matter, etc. [7j. Moreover, 
they are seriously constrained by the tight limits on the photon fraction from the latest Auger 
measurements |12]. In the following, we concentrate on the bottom-up, astrophysical scenarios. 

A spectral steepening analogous to that for protons can also result if the highest energy 
UHECRs are composed mainly of heavy nuclei such as iron. For nuclei, the dominant energy 
loss process above 10^^ eV during intergalactic propagation is photodisintegration and pair 
production interactions with photons of the far-infrared background (FIRB) and the CMB [13j. 
Based on recent determinations of the FIRB, the energy loss length for iron nuclei at 10^'^ eV is 
^ 300 Mpc, somewhat larger than that for protons [T3|. Observationally, the nuclear composition 
of UHECRs, particularly at the highest energies, is quite uncertain. Fluorescence measurements 
of shower elongation rates by the HiRes experiment have indicated that the composition changes 
from heavy-dominated below 10^^ eV (presumably of Galactic origin) to light-dominated above 
this energy [15]. However, the latest results from Auger, based on hybrid events with higher 
precision and considerably higher statistics, instead reveal an intermediate composition at all 
energies up to 4 x 10^^ eV, with hints of a trend toward a heavier composition at the highest 
energies [IB]. A caveat is that systematic uncertainties in interaction models and atmospheric 
optical attenuation are still quite large [TTj. These issues are also critical for the interpretation 
of the ankle feature in the spectrum at ~ 3 X 10^^ eV [ig [I9]. The Telescope Array should 
provide key complementary information in this regard, since its scintillator surface detectors 
are relatively more sensitive to the electromagnetic component of the air shower than the muon 
component, in comparison with the water Cherenkov tanks of Auger. 

Another unavoidable aspect is deflection by Galactic and extragalactic magnetic fields 
(EGMF), which affect the arrival directions of UHECRs and can also act to lengthen their 
effective propagation distance. The strength and distribution of magnetic fields in the 
intergalactic medium as well as at high latitudes in the Galaxy are very poorly known, both 
observationally and theoretically. Faraday rotation measurements of distant radio sources 
give only upper limits in the nanogauss range for intergalactic fields on average, subject to 
assumptions on the field coherence scale and ionized gas distribution L20J . Realistically, whatever 
their origin, EGMF can be expected to have some correlation with the distribution of large scale 
structure, and various EGMF models that account for this have been introduced to assess their 
effects on UHECR propagation [21J. The effect of Galactic magnetic fields (GMF) is also model- 
dependent [^1^. 

The Auger collaboration has recently announced statistically significant evidence for UHECR 
anisotropy [H] . They find a ~ 3(T correlation between the arrival directions of 20 UHECR events 
above 5.6 x 10^^ eV, and the positions of active galactic nuclei (AGNs) in the Veron-Cetty and 
Veron catalog with distances <75 Mpc, to within ~ 3 degrees. Some robust implications of this 
result are that: 1) the UHECR sky is not isotropic, 2) the sources of UHECRs are extragalactic 
and trace large-scale structure on scales of order 100 Mpc, ruling out Galactic sources, and 3) 
deflections of UHECRs in the GMF and EMGF are not so severe as to isotropize the arrival 
directions, so that charged particle astronomy is possible. 

Not much more can be said with confidence, however. A correlation with the global 
distribution of AGNs does not necessarily mean that they are the true sources; any other type of 
source that follows the distribution of large-scale structure, such as gamma-ray bursts (GRBs) or 
clusters of galaxies, may also be feasible. Inferences on the source density or UHECR composition 
from the current anisotropy results are dependent on assumptions regarding the EGMF and 



GMF (see below). In addition, the Veron-Cetty catalog is known to be highly non- uniform 
and incomplete, and further studies with higher statistics utilizing more suitable catalogs are 
warranted (see e.g. [25] for recent studies in this direction). Since the northern hemisphere 
should be less affected by the GMF, data from the Telescope Array will be crucial, and all-sky 
observations by JEM-EUSO even better. We note that the AGN correlation results have not 
been confirmed by the northern HiRes, although their statistics is smaller |26] . At present, the 
true identity of the UHECR sources is still wide open. 

3. Acceleration and energetics in astrophysical sources of UHECRs 

A minimum requirement for astrophysical sources of UHECRs is the ability to magnetically 
confine particles of the requisite energies. For particles with energy E and charge Z, this implies 
the Hillas condition (i?/pc)(S/lG) ;^ {E/10'^'^eY)/Z between the system's size R and magnetic 
field B [27]. Note that this can also be rewritten as a lower limit on the source power, assuming 
that the magnetic field carries a fixed fraction of the outfiowing energy density \28\ |29] . Only a 
few types of objects are known to meet this criterion, among them the jets of radio-loud AGNs, 
GRBs, and clusters of galaxies. Below we focus on these three as representative types of potential 
UHECR sources (see [27] for other possibilities). The actual maximum energy attainable under 
different circumstances must be evaluated by comparing the timescales for particle acceleration, 
often that for the first-order Fermi mechanism in shocks, against the timescales for limiting 
processes such as source lifetime, particle escape, adiabatic or radiative energy loss, etc. 

Equally important is the available energy budget. FigH] shows estimates of the kinetic energy 
output averaged over the universe as a function of redshift z due to AGN jets, GRB explosions 
and accretion onto clusters, which should be roughly proportional to their cosmic ray output. 
The plotted quantity is differential per unit z, dE\f_i^/dz = (dt/dz) J L{dn/dL)dL, where L 
is the kinetic luminosity per object and dn/dL is the z-dependent luminosity function, with 
cosmological parameters h=0.7, r2m=0.3 and Q\=0.7. For AGN jets, we have made use of the 
observed radio luminosity function along with the observed correlation between the radio and 
jet kinetic luminosities of radio galaxies [30J. GRBs were assumed to occur each with kinetic 
energy -Bgrb = 10^^ erg at a rate that follows the star formation history and matches the 
log iV-log S distribution observed by BATSE ^31j , the estimate being roughly independent of 
the beaming factor. Evaluations based on the more recent, post-Swift z-distribution [32] are 
also shown. The different curves for AGNs and GRBs correspond to different evolutionary 
assumptions at the highest z, with only small differences at low z. For clusters with mass 
M, the rate of gas kinetic energy dissipation through accretion shocks can be estimated as 
Li,cc ^ 9 X lO^^(M/lO^^M0)5/3 erg/s |331|3l]. This can be combined with the Press-Schechter 
mass function to evaluate dE-^i^^/dz for clusters of different M. Note that due to the hierarchical 
nature of structure formation together with the nonlinear nature of gravity, the maximum is 
reached at z = 0. 

The results at low z can be compared with the observed energy density of UHECRs, 
~ 10~^^erg cm~'^ ~ 3 x 10^'^erg Mpc~^ above 10^^ eV. Considering further factors for energy 
loss during propagation, it is apparent that whereas AGN jets and cluster accretion shocks have 
reasonable margins to accommodate the energetics of UHECRs, GRBs, with a substantially 
smaller energy budget, require a high efficiency of energy conversion into UHECRs. 

4. Active galactic nuclei 

AGNs are luminous active objects in the centers of some galaxies that are powered by accretion 
of matter onto supermassive black holes [33]. Roughly 10 % of them are known to be radio- loud, 
ejecting powerful, relativistic jets from the nucleus out into its surroundings, but the majority 
of AGNs are of the radio-quiet type with no strong jets. In contrast to the former, the latter 
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Figure 1. Energy budget of 
candidate UHECR sources: AGN 
jets (dotted), GRBs (dashed; thin 
for BATSE, thick for Swift) and 
cluster accretion shocks (sohd), the 
latter separately for each logM as 
labelled. See text for more details. 
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are considered unlikely to be UHECR sources due to the severe radiative losses expected near 
the nucleus [36j, and the absence of nonthermal emission components. 




Radio-loud AGNs themselves come in two basic types [37]. In FR II objects with radio (or jet 
kinetic) power above a critical value, the jets remain relativistic out to large scales and end in 
luminous hot spots, termination shocks where the jets are strongly decelerated upon impacting 
the intergalactic medium. If such jets are oriented close to the line of sight, they would appear as 
GeV-bright, flat spectrum radio quasars. Contrastingly, the jets of FR I objects with subcritical 
power are initially relativistic but is believed to decelerate gradually as they propagate out of 
their host galaxies. Their on-axis counterparts are the TeV-bright, BL Lac objects. 

Different locations along AGN jets can be UHECR production sites, since one expects 
B oc under the naive assumption that the ratio of magnetic to kinetic energy is constant. 
For both FR II and FR I objects, a candidate is the inner jet region with R ~ 10^^-10^^ cm and 
B ~ 0.1-1 G, known through observations of blazars to be a site of intense particle acceleration, 
perhaps due to internal shocks. Estimates assuming Bohm-type acceleration (i.e. acceleration 
time comparable to gyration time in fully turbulent magnetic fields) show that the maximum 
proton energy should be limited by photopion interactions with low frequency internal radiation 
to somewhat below lO^'' eV p8]. However, plausible models of broadband blazar emission based 
on electron acceleration generally indicate that the plasma conditions in inner jets are far from 
the Bohm limit [39]. Furthermore, conversion to neutrons may be necessary for the particles 
to escape the jet without suffering adiabatic expansion losses and contribute to UHECRs by 
decaying back to protons outside. This may entail rather finely-tuned conditions, as photopion 
interactions must be efficient to yield enough neutrons, yet cannot be excessively so to avoid 
adverse energy losses. A less problematic site may be the hot spots of FR II radio galaxies with 
R ~ lO'^^ cm and i? ~ 1 mG, where the maximum energy may exceed 10^*^ eV, limited by escape 
|4Uj . However, note that these particles must further traverse the extensive cocoon of shocked, 
magnetized jet material in order to completely escape the system, an issue that has not been 
examined in detail. A further possibility is acceleration by the bow shocks being driven into the 
external gas by the expansion of the cocoon, if the ambient magnetic field can be sufficiently 
amplified [28j. In view of the Auger composition results [IH], another interesting question is 
whether heavy nuclei such as iron can be accelerated and survive photodisintegration in these 
environments [4lj . 

In order to verify an AGN origin, detailed analysis of observed UHECR arrival directions and 
cross correlations with known source positions is undoubtedly essential, and the recent Auger 




results are an important step in this direction. However, it is not a trivial task, given the 
uncertainties in the intervening GMF and EGMF and the UHECR composition (Sec. [2]). In 
addition to deviations in the arrival directions, deflections in the EGMF during propagation also 
entail significant delays in the arrival times of UHECRs, which can be of order 10'' yr or more for 
distances 100 Mpc, depending on the properties of the EGMF [42j. This is in fact comparable 
to or longer than the typical lifetime or duration of AGN activity, particularly for radio-loud 
AGNs [43j, so that they may effectively behave as transient sources, in a way similar to GRBs. 
Before the UHECRs emitted by an AGN reaches the observer, the object may significantly evolve 
in its radio power, or possibly shut off its activity altogether. Another implication is that an 
effectively larger number of sources contribute to small scale clustering, underscoring the need 
for high event statistics as might be provided by JEM-EUSO or Auger North+South. 

Because of these concerns, it is also highly desirable to have some means to pinpoint 
individual sources through characteristic, UHECR-induced signatures of secondary neutral 
radiation. Neutrinos are ideal in providing an unambiguous earmark of high energy hadrons 
[44j . However, even with km^ detector facilities such as IceCube or KMSNeT, it may not be 
easy to study individual AGNs in sufficient detail. Thus, distinctive electromagnetic signals 
will also be extremely valuable. For radio galaxy hot spots, synchrotron emission from UHE 
protons can produce nonthermal X-rays that could be distinguished from other processes such as 
electron inverse Compton through multiwavelength observations [45j . If the medium surrounding 
the source is sufficiently magnetized, UHE protons propagating in the source's vicinity can lead 
to diffuse gamma-ray emission from photomeson-triggered cascades that may be detectable by 
current and upcoming instruments |46j . 

5. Gamma-ray bursts 

GRBs are explosive phenomena that are believed to arise in a fraction of stellar collapse events, 
characterized by highly variable, ultrarelativistic outflows lasting for ~ 10 — 100 s [4T| . As with 
AGN jets, these outflows can be host to different sites for UHECR acceleration. Potential locales 
include internal shocks, external reverse shocks, and external forward shocks, believed to be the 
emission sites of the prompt X-rays and gamma-rays, optical flash and radio flare, and the radio 
to X-ray afterglow, respectively [48j. The external forward shock may possibly be disfavored due 
to the ultrarelativistic velocity and weak upstream magnetic field |;49j, but this is controversial 
|50j . For the mildly relativistic internal and external reverse shocks, a different problem is that 
for the particles to escape the acceleration site without significant losses, neutron conversion 
may be required, as with AGN inner jet regions. 

A crucial distinguishing feature of UHECRs from GRBs is the narrow dispersion in CR energy 
expected for individual sources, due to the time delay during propagation in EGMF [32]. Clear 
detections of this effect demand substantial event statistics that may only be achievable with 
future facilities. Even then, it would only demonstrate that UHECRs come from bursting 
sources, and their GRB origin will remain ambiguous. Thus, multi-messenger photon and 
neutrino signatures are particularly essential in the case of GRBs. 

Identification of high energy neutrino signals from GRBs are facilitated through time 
coincidence, but studying individual bursts in detail may be difficult. Thus, photon signatures 
of UHECR production will also be crucial [51] ■ Detailed studies of this issue were conducted 
for the GeV-TeV emission from internal shocks, utilizing a comprehensive Monte Carlo code 
that includes a wide variety of physical processes related to high energy protons and electrons 
\52\ I53j . Besides inverse Compton emission from primary electrons, interesting proton- 
induced components can become clearly visible, such as synchrotron emission from high energy 
protons and muons, as well as emission from secondary pair cascades injected by photomeson 
interactions, serving as unique signatures of ultra-high-energy protons. 

Particularly interesting are proton-dominated GRBs that contain a significantly larger 



amount of energy in accelerated protons compared to accelerated electrons |53j. Such conditions 
are motivated not only by the physics of particle acceleration in collisionless shocks, but also by 
the latest measurements of the GRB z-distribution. Post-Swift estimates of the local GRB rate 
\32\ [5l] indicate that the proton energy content per burst may need to be substantially higher 
than previously thought in order for them to remain viable as UHECR sources. Unique UHECR- 
induced GeV-TeV components are then expected, such as GeV peaks, UV- X-ray excesses and 
luminous TeV bumps (FigO]). 
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Figure 2. Prompt emission spectra from internal shocks in proton-dominated GRBs (see [53] 
for more details). (Left) For F = 300 and different proton/electron ratios. Solid curves are 
the total emission, and dot-dashed curves denote separately electron synchrotron (eSY) and 
electron inverse Compton (elC) components without 77 absorption effects. Dashed curves are 
primary electron components only. (Right) Similar to left panel, except for F = 1000 and 
different magnetic fields. Dot-dashed curves denote separately proton synchrotron (pSY) and 
muon synchrotron (^SY) components. 

Note that the conditions favorable for GeV- TeV emission may also imply efficient UHECR 
acceleration and escape, but not necessarily strong neutrino emission |52] . The observational 
prospects are promising for GLAST, atmospheric Cherenkov telescopes such as MAGIC (II), 
H.E.S.S. (II), VERITAS and CANGAROO III, as well as wide-field surface detector facilities, 
which should test these expectations and provide important information on the physical 
conditions in GRB outflows. 

A different class of low luminosity GRBs have been proposed to be more important than 
high luminosity GRBs for the integrated energy budget, and hence potentially more relevant as 
sources of UHECRs [55J. As with AGN jets, the acceleration and survival of heavy nuclei in in 
the environments of both high and luminosity GRBs is an important question [56j . 

6. Cluster accretion shocks 

Clusters of galaxies are the largest gravitationally bound systems in the universe, composed 
of stars, hot gas and dark matter. According to the currently favored picture of hierarchical 
structure formation in the CDM cosmology, they are the latest objects to form, and should be 
surrounded by strong accretion shocks as a consequence of continuing infall of dark matter and 
baryonic gas [58]. Such shocks should be interesting sites of particle acceleration and have been 
proposed as sources of UHECRs [28] . However, estimates of the maximum energy -Emax for 
protons seem to fall short of lO^'^ eV by 1-2 orders of magnitude (SHI [Ml- 



We have shown that invoking UHECR nuclei in cluster accretion shocks may offer a possible 
solution [57]. Heavy nuclei with higher Z have correspondingly shorter acceleration time so that 
Fe may be accelerated up to 10^" eV in the same shock conditions, notwithstanding energy losses 
by photodisintegration with the FIRB and CMB. With reasonable assumptions regarding the 
number density and CR power of the sources, detailed propagation calculations demonstrate that 
the model compares favorably with the existing observations, as long as the source composition 
of Fe is similar to that of Galactic CRs (FigJS]). See also 1^0] for related studies on UHECR 
nuclei propagation. 
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Figure 3. Observed UHECR spectrum (left) and mean mass composition (right) versus energy 
E (1 EeV = 10^^ eV) from cluster accretion shocks for a = 1.7 and (3 = 0.5 [57|. Compared are 
current data for HiRes [61J (bars) and Auger [TU] (asterisks). The histograms are the average 
over different model realizations for the cases with (thick) and without (thin) EGMF, and the 
thin curves outline the median deviations due to cosmic variance for the former case only. The 
straight line in the top panel denotes the injection spectrum. 



The mass composition at ^ 3 x 10 eV should be predominantly light and consistent with 
HiRes reports [.15] . The rapid increase of the average mass at higher energies is a clear prediction 
of the scenario, in line with the latest Auger elongation measurements [16], and to be tested 
further by new generation facilities including Telescope Array and JEM-EUSO. 

Despite the relative rarity of massive clusters in the local universe, strong deflections of 
the highly charged nuclei in EGMF allow consistency with the global isotropy. On the other 
hand, with a sufficient number of accumulated events, anisotropics toward a small number of 
individual sources should appear, although this expectation is subject to uncertainties in the 
EGMF and GMF. Whether it is consistent with the latest Auger anisotropy results calls for 
further investigation. It could possibly be the case, including the observed deficit of events from 
the Virgo direction, if the GMF is of the plausible, bisymmetric spiral type, in which deflection 
angles can be minimal in selected regions of the sky around the supergalactic plane, particularly 
around Centaurus [23] (see also [62] for related discussions). Note that the Centaurus cluster 
at distance ~50 Mpc is one of most prominent clusters in the local universe, somewhat more 
massive than Virgo [63j . 

We mention that the escape of CRs from clusters can be mediated by diffusion in directions 
away from the filaments, or perhaps by advection in partial outflows during merging events. 
Energy-dependent escape upstream of the shock may also be possible. An aspect of this scenario 
that warrants further study is the spectral domain < 10^^ eV and the implications for the 
Galactic-extragalactic transition region. Alternatively, provided that EGMF horizon effects are 
effective ^ 10^^ eV [6l], cluster accretion shocks may give an appreciable contribution to CRs 



between the second knee and the ankle, in which case they may also be interesting as high energy 
neutrino sources [65] . 

Finally, we may look forward to very unique signatures in X-rays and gamma-rays. Protons 
accelerated to 10^^-10^^ eV in cluster accretion shocks should efficiently channel energy into 
pairs of energy 10^^-10^^ eV through interactions with the CMB, which then emit synchrotron 
radiation peaking in hard X-rays and inverse Compton radiation in TeV gamma-rays [34J. The 
detection prospects are promising for Cherenkov telescopes such as HESS, and hard X-ray 
observatories such as Suzaku and the future NeXT mission. Photopair production by nuclei 
may also be efficient and induce further interesting signals that are worth investigating. 

7. Conclusion and outlook 

Besides AGNs, GRBs and clusters, other proposals for the astrophysical origin of UHECRs not 
covered in this brief review include starburst galaxies [66] , extragalactic magnetars |67j , dormant 
black holes [68], etc. (see [27]). The true source of UHECRs could be one or more of these; 
alternatively, we should not discount the possibility that it is actually none of them, and the 
truth something totally unexpected. 

In this regard, we may recall the history of GRB research. Before the launch of BATSE 
in 1991, the sample of GRBs numbered ^ 200, with no clear anisotropy observed in the sky 
distribution. However, most believed, based on theoretical plausibility, that they were generated 
by Galactic neutron stars, while a very small minority advocated neutron star mergers at 
cosmological distances. After several years and a sample of more than 3000 bursts, the currently 
most favored scenario for the progenitors of long GRBs is related to the collapse of massive stars 
[69j , a concept that didn't exist before BATSE. The lesson is that we should be cautious about 
jumping to conclusions on the origin of UHECRs based on just 20 observed events and our 
limited theoretical understanding of cosmic objects. 

Regardless of theoretical prejudice, the quest for the solution of the UHECR mystery is 
progressing in earnest with ongoing measurements by Auger and Telescope Array. An order of 
magnitude greater statistics for the highest energy events expected from JEM-EUSO or Auger 
North-I-South, and the combined effort of upcoming neutrino. X-ray and gamma-ray observations 
is expected to lead us toward the true answer in the near future. 
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